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We propose a method for Monte Carlo simulations of systems with a complex action. The method has the 
advantages of being in principle applicable to any such system and provides a solution to the overlap problem. 
In some cases, like in the IKKT matrix model, a finite size scaling extrapolation can provide results for systems 
whose size would make it prohibitive to simulate directly. 



1. INTRODUCTION 

There exist many interesting systems in high 
energy physics whose action contains an imag- 
inary part, such as QCD at finite baryon den- 
sity, Chern-Simons theories, systems with topo- 
logical terms (like the #-term in QCD) and sys- 
tems with chiral fermions. This imposes a severe 
technical problem in the simulations, requiring an 
exponentially large amount of data for statisti- 
cally significant measurements as the system size 
is increased or the critical point is approached. 
Furthermore, the overlap problem appears when 
standard reweighting techniques are applied in 
such systems and it becomes exponentially hard 
with system size to visit the relevant part of the 
configuration space. In it was proposed to take 
advantage of a factorization property of the dis- 
tribution functions of the observables one is in- 
terested to measure. This approach can in prin- 
ciple be applied to any system and it eliminates 
the overlap problem completely. In some cases 
it is possible to use finite size scaling to extrap- 
olate successfully to large system sizes where it 
would have been impossible to measure oscillat- 
ing factors directly. The method has been applied 
successfully in matrix models of non perturbative 
string theory (IKKT) [Jj, random matrix theory 



of finite density QCD (RMT) Z as well as the 
2d CP 3 model, the Id antiferromagnetic model 
with imaginary B and the 2d compact U(l) with 
topological charge [3]. 

In this paper we present our results for 
IKKT and RMT. In the first case we study the 
space-time dimensionality hoping to dynamically 
recover our 4d space-time and in the second to 
test the factorization method against known an- 
alytical results. In all cases we will be deal- 
ing with a system defined by a partition func- 
tion Z = f dAe~ s ° e lF and the corresponding 
phase quenched model Zq — J dAe~ s ° where 
S = Sq—iT is the action of the system with its real 
and imaginary parts. A represents collectively the 
degrees of freedom of the model and in our case it 
corresponds to a set of N x N matrices. In case 
we are interested in measuring some observable 
O, we consider the distribution functions po{x) — 
(6(x-0)) and p% ] '(a?) = (5(x-O)) , where (.. .)„ 
refers to Zq. Then we define the fiducial sys- 
tem Zq.x = J dAe~ So 6(x — O), the weight fac- 
tor wq(x) = (e tT )o,x and the distribution po(x) 
factorizes po{ x ) = ^jPo( x ) w o( x ) where C = 
<e* r } . Then (O) = £ dx x p% ] (x) w a (x) . 
The (5-function constraint is implemented in our 
simulations by considering the system Zq v = 



2 



JdAe- So e y <°) where V(z) = ±7(2 -£) 2 and 7 ,£ 
are parameters which control the constraining of 
the simulation. The results are insensitive to the 
choice of 7 as long as it is large enough. Then we 
have that wq(x — (0)iy) = (e tT )i.v- The distri- 
bution of O in Ziy has a peak x and the quan- 
tity V'{x) is the value of f$\x) = £\np^\x) 

at x = x. The function /^(rc) can be obtained 
by integrating an analytic function to which we 
fit the /q (x) data points. 

By applying this method we force the system 
to sample configurations which give the essential 
contributions to (O) , something that would be in- 
creasingly difficult with system size in the phase 
quenched model, eliminating this way the over- 
lap problem. This already allows us get close to 
the thermodynamic limit with modest computer 
resources. Furthermore we obtain direct knowl- 
edge of wq{x) and po(%) which allows us to un- 
derstand the effect of T. This is important for 
understanding the properties of the system when 
r plays a crucial role. Using the generic scaling 
properties of the weight factor Wi(x), one may 
extrapolate the results obtained by direct Monte 
Carlo evaluations to larger system size. Such an 
extrapolation is expected to be particularly use- 
ful in cases where the distribution function turns 
out to be positive definite. In those cases we can 
actually even avoid using the reweighting formula 
by reducing the question of obtaining the expec- 
tation value to that of finding the minimum of 
the free energy, which is (minus) the log of the 
distribution function. Here, the error in obtain- 
ing the scaling function propagates to the final 
result without significant magnifications. There- 
fore, the extrapolation can be a powerful tool to 
probe the thermodynamic limit from the accessi- 
ble system size. 

2. NON PERTURBATIVE STRING 
THEORY 

The so-called IKKT matrix model 0] has 
been conjectured to be a definition for non- 
perturbative string theory. A particularly inter- 
esting feature of the model is that the eigenval- 
ues of the bosonic matrices generate space- 



time dynamically. In our case we work with a 
reduced model where /1 = 1, . . . , 6. The observ- 
ables that we study are the normalized eigenval- 
ues Ai/(Ai) > ... > A 6 /(A 6 ) of the space-time 
"moment of inertia" — -^TrA^A^. An in- 
teresting scenario to investigate would be that 
the 0(6) symmetry is spontaneously broken with 
some of the eigenvalues (possibly 4) grow large 
and the rest remain small, providing a mecha- 
nism for dynamical compactification of extra di- 
mensions. 
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Fig. 1: logp 5 (a;) for N=128, n = 16. 

The IKKT model is a good example where fi- 
nite size scaling for the oscillating factor Wi(x) 
works well (the index i corresponds to Ai/(A;}o). 
For i > 1 we find fast convergence to a scal- 
ing function <f>i(x), where Wi(x) = exp{N 2 $i(x)} 
PP. We compute &i(x) for small matrix size n < 
20 and the phase quenched distribution function 
pf\x) for larger size N. Then the factorization 
formula can be used in order to compute pi{x). 
Note that in the computation of (Aj) the errors 
do not propagate exponentially with system size 
since (Aj) can be determined from the minimum 
of the "free energy" Fi(x) = — -A^ log pi(x). In 
Fig. 1 we show our results for ps(x) for n = 16 
and N = 128. A double peak structure of ps(x) 
is evident and it is expected that the peak for 
small x will increase with system size £Q. One 
hopes that this peak will be dominant for i = 5, 6 
and that the large x peak will be dominant for 
i = 4, 3, 2, 1 and this will realize the SSB sce- 
nario. Note the heavy suppression of the x = 1 
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region caused by T which is the peak of the phase 
quenched model distribution of A;. 
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values of N < 48 to obtain the thermodynamic 
limit. Unfortunately, the function wr(x) is not 
positive definite and the important contributions 
come from the region where it changes sign. As 
expected, we find that finite size scaling does not 
work as well as in the case of the IKKT model 
(although we obtain agreement up to order of 
magnitude for the values of N < 96 that we 
explored). We also find it very difficult to ex- 
plore the crossover region near the phase transi- 
tion point fi c = 0.527 . . . for N > 8 since |tUi(x)| 
becomes very small. Since RMT is a schematic 
model of finite density QCD, we expect that the 
factorization method will be useful to explore the 
phase diagram of QCD. 
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Fig. 2: Pr(x) and pi{x) for /i = 0.2. 



Fig. 3: (u) for N = 8. 



3. RMT OF FINITE DENSITY QCD 

We consider RMT with one quark flavour and 
zero quark mass [Jy. The model is chosen in or- 
der to study the correctness and effectiveness of 
the factorization method, since one can compare 
results with known analytical solutions even for fi- 
nite N. The observable we measure is the "quark 
number density" v as a function of the chemi- 
cal potential //, and we consider the distribution 
functions pi(x), where i = R,I corresponds to 
the real and imaginary parts of v respectively. 
Notice that the effect of T is dramatic, causing a 
discontinuous transition in v. In [2] our results 
nicely reproduce the exact results known for fi- 
nite N and we are able to achieve large enough 
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